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ABSTRACT 


A  unified  theory  on  a  certain  type  of  Inequalities  for 
abstract  linear  operators  and,  in  particular,  for  ordinary 
linear  differential  operators  of  the  first,  second,  and  fourth 
order  is  developed.  The  statements  which  involve  these 
inequalities  lead  to  a  principle  of  error  estimation.  With  a 
programmed  procedure,  approximate  solutions  and  corresponding 
error  bounds  have  been  calculated  for  some  examples  of  the 
second  and  fourth  order. 


1.  INTRODUCTION 


This  paper  is  mainly  concerned  with  linear  ordinary  differential 
operators  L  such  that 

L[u](x)  21  0  (0  <  x  <  1)  implies  u(x)  _>  0  (0  <.  x  <_  1) 

for  all  functions  u(x)  which  possess  derivatives  of  sufficient  high 
order  and  satisfy  certain  equations ,  or  inequalities ,  at  the  boundary 
points  x  ■  0  and  x  ■  1.  In  particular ,  we  are  interested  in  how 
this  property  can  be  applied  in  order  to  obtain  error  bounds  for  an 
approximate  solution  of  a  boundary  value  problem  which  involves  the 
operator  L. 

In  Section  2,  an  abstract  theory  on  inverse-positive  linear  operators 
M  is  developed.  This  theory  is  then  applied  to  differential  equations 
of  the  first  and  second  order  (Section  3),  to  systems  of  such  equations 
(Section  4),  and  to  some  differential  equations  of  the  fourth  order 
(Section  5).  For  all  of  these  problems,  we  obtain  in  this  way  principles 
for  error  estimation.  In  Section  6,  we  finally  give  some  numerical  results 
which  have  been  gained  by  a  programmed  procedure  based  on  the  error 
estimation  principle  just  mentioned. 

Many  results  concerning  inequalities  for  ordinary  linear  differential 
operators  can  be  found  in  the  literature.  In  order  to  prove  these 
results,  several  different  methods  have  been  used.  This  paper  shall  show 
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that  many  of  those  results,  and  new  results  also,  can  be  derived  In 
a  rather  simple  way  from  simple  abstract  theorems.  In  particular, 
this  is  true  for  results  which  are  useful  for  numerical  calculations. 
Some  of  the  theorems  in  this  paper  are  already  known,  for  example,  the 
basic  Theorem  1  in  Section  2,  and  its  applications  to  differential 
operators  of  first  and  second  order,  and  also  the  remarks  concerning 
iterative  procedures  in  Section  2.6.  Most  of  the  other  results  seem 
to  be  new. 

The  results  in  this  paper  can  be  generalized.  For  example,  the 
abstract  theorems  in  Section  2,  as  well  as  the  numerical  procedure  in 
Section  6,  can  be  applied  to  linear  partial  differential  equations  also. 
Moreover,  Theorem  1  in  Section  2  is  a  special  case  of  a  theorem  on 
nonlinear  inverse-monotonic  operators  M.  This  more  general  theorem 
can  be  applied  to  nonlinear  ordinary  and  partial  differential  equations, 
and  a  similar  procedure  as  that  in  Section  6  can  be  used  to  obtain  error 
bounds  for  certain  types  of  such  nonlinear  equations.  While  a  survey 
on  such  applications  to  nonlinear  equations  will  be  given  elsewhere  [8] , 
we  restrict  ourselves  here  to  linear  operators  and  ordinary  differential 
equations  in  order  to  be  able  to  give  more  details. 
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2.  INVERSE-POSITIVE  OPERATORS 

Let  R  ■  {u,v,...}  and  S  ■  {U,V,...}  be  linear  partially 
ordered  spaces,  and  let  M  be  a  linear  operator  which  naps  R  into  S. 
Under  what  conditions  is  it  true  that 

Mu  >_  0  implies  u  >_  0 

for  u  e  R? 

A  linear  operator  with  this  property  will  be  called  inverse - 
positive ,  because  the  above  implication  holds  if  any  only  if 

M  has  a  positive  inverse , 

i.e.,  M  ^  exists,  and  U  0  implies  M  >.  0  for  U  in  the 
domain  of  M 

2.1  SUFFICIENT  CONDITIONS 

Let  R  be  Archimedian,  i.e. 

nu  v  (n  *  1,2 , . . . )  implies  u  0 

if  u,v  are  arbitrary  fixed  elements  in  R.  We  define  a  second  order 
relation  u  v  in  R  (  u  strongly  smaller  than  v  )  by: 

u  yo  (  u  strongly  positive)iff  u  has  the  following 

PROPERTY  o:  For  each  v  e  R  there  exists  a  natural  number  n 
such  that  v  <_  nu. 

Moreover,  we  require  that  there  is  defined  a  second  order  relation 
V  in  S,  also.  This  relation  may  be  defined  by  Property  o,  or  may 
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not,  but  it  shall  satisfy  the  following  conditions: 

U>-0,  V  >_  0  imply  U  +  vX'O;  U>»0,  A>0  imply  AU>-0; 

U  ^0  implies  U  >_  0. 

These  assumptions  shall  hold  throughout  this  Section  2. 

THEOREM  1:  Let  the  following  assumptions  be  satisfied. 

Assumption  I:  For  w  e  R, 

w  ^  0,  Mw  ^*0  imply  w  ^*0. 

Assumption  II:  There  exists  z  e  R  such  that 

z  >  0,  Mz  ^-0.  (2.1) 

Then  for  u  e  R, 

Mu  >_  0  implies  u  >_  0. 

PROOF  (see  Figure  1):  Suppose  that  Mu  ^  0,  but  not  u  0. 
Then,  because  R  is  Archimedian,  there  exists  a  smallest  number 
AQ  >  0  such  that  u  +  A^  z  >_  0.  The  element  w  *  u  +  AQz  satisfies 
v  >_  0  but  not  w  ^»0.  On  the  other  hand,  Assumption  II  yields 
Mw  *  Mu  +  AqMz  ^"0,  and  therefore,  w  ^0  because  of  Assumption  I. 
This  contradiction  shows  that  the  Theorem  is  true.  (Notice,  that 
z  >  0  because  of  Assumption  I.  For  a  more  detailed  proof,  see  [7].) 


Figure  1:  Illustration  of  the  proof  with  R  ■  S  - 

12  1  2 
u  *  (u  ,u  )  >.  0  iff  u  ^  0,  u  >  0;  u>-0  iff 


and  the  same  order  relations  in  S. 
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SPECIAL  CASE  1:  If  M  ■  A-B  with  a  linear  operator  A 
satisfying  Assumption  I,  and  a  positive  linear  operator  B,  then 
M  satisfies  Assumption  I.  In  particular,  A  «  I  satisfies 
Assumption  I. 

2.2  NECESSARY  CONDITIONS 

We  will  prove  that  the  sufficient  conditions  of  Theorem  1  are 
necessary  "in  general". 

THEOREM  2:  Suppose  that  M  is  inverse  -positive .  Then  the  follow¬ 
ing  statements  hold. 

1.  The  operator  M  satisfies  Assumption  I  if  the  second  order 
r  elation  in  S  has  Property  o- 

2.  The  operator  M  satisfies  Assumption  II  if  its  range  contains 
at  least  one  strongly  positive  element. 

PROOF:  1)  Suppose  that  w  >_  0,  Mw  ^  0  for  some  we  R  and  let 

v  be  any  element  in  R.  Then,  if  the  second  order  relation  in  S  has 
Property  o,  there  exists  a  number  n  satisfying  Mv  nMw.  This 
inequality  yields  v  <_  nw  for  inverse-positive  M.  Because  there  exists 
such  a  number  n  for  every  v  e  R,  we  have  w  ^*0  by  definition. 

2)  Let  rcS  be  any  strongly  positive  element  in  the  range  of  M  and 
Mz  *  r.  Then  Mz  >_  0  and  therefore  z  ^  0,  because  M  is  inverse¬ 
positive.  Thus,  z  satisfies  Assumption  II. 
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2.3  CONNECTED  CLASSES  OF  OPERATORS 

For  many  important  classes  of  operators,  it  is  easy  to  prove 
that  Assumption  I  is  satisfied  (see,  for  example,  Sections  3  and  4). 

If  M  is  inverse-positive,  the  element  z  occurring  in  Assumption  II 
may  be  constructed  by  solving  the  equation  Mz  *  r  for  any  r  0 
(see  the  second  part  of  the  proof  of  Theorem  2).  Considering,  however, 
certain  "connected"  classes  of  linear  operators  M,  instead  of  a  single 
operator,  one  need  prove  Assumption  II  for  only  one  operator  of  such  a 
class. 

Let  9)?  be  a  set  of  linear  operators  M  which  map  R  into  S 
and  satisfy  Assumption  I. 

A  one  parameter  family  z(t)  e  R  (0  <_  t  1)  will  be  called 
continuous  at  t^  e  (0,1)  if  for  each  natural  n  there  exists  a  positive 
number  6(n)  such  that 

-~e  <_  z ( t )  -  z(tQ)  <_~e  for  |t-tQ|  1  6(n),  (2.2) 

where  e  0  denotes  an  arbitrary  element  in  R  (which  may  depend 

on  z  and  t  ) . 

0 

LEMMA:  Suppose  that  for  each  pair  of  operators  M^,  M  e  9)1 
there  exists  a  family  of  operators  M(t)  e  2)?  and  a  continuous  far'd  ly 
z(t)  e  R  (0  t  <_  1) ,  such  that 

M(0)  -  M0,  M(l)  -  Mr  M(t)z(t)>-0  (0  <  t  <  1). 

lety  moreovery  at  least  one  operator  M  e  9)?  be  inverse-positive. 
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Then  all  operators  in  SR  are  inverse-positive. 

PROOF:  Suppose  c  SR  is  not  inverse-positive.  Define 
MQ  *  M  and  let  M(t)fz(t)  have  the  properties  mentioned  above. 

The  inequality  Mqz(0)  ^ 0  implies  z(0)  ^0  because  MQ  is 
inverse-positive.  The  two  last  inequalities  yield  z(0)  ^-0  because 
of  Assumption  I.  Thus,  there  exists  a  smallest  positive  number 
tQ  <_  1  such  that  z(t)  ^0  for  0  t  <  tQ.  For  this  number  tQ, 
the  right  hand  inequality  in  (2.2)  yields 

n(-z(tQ))  <_  e  -  nz(t)  <_  e  if  0  <  tQ  -  6  <  t  <  tQ. 

Therefore,  we  have  -z(tQ)  0  because  R  is  Archimedian. 

If  tQ  ■  1,  the  element  z(l)  satisfies  (2.1)  for  M  ■  and, 
therefore,  is  inverse-positive  according  to  Theorem  1.  The 

inequality  t^  <1,  however,  contradicts  the  definition  of  t^. 

Applying  Assumption  I  to  M(tQ),  we  get  z(tQ)  Therefore,  there 

exists  a  natural  number  m  with  the  property:  e  <^mz(t0).  For  n  >  m, 
the  left  hand  inequality  in  (2.2)  yields 

0«  (1  -  J&)z(t  )  <_  z(t)  if  |t  -  tQ|  small  enough. 

That  means,  z(t)  ^0  is  true  for  a  larser  interval  than  [0,tQ). 

One  need  not  always  actually  construct  the  elements  z(t).  In  order  to 
prove  their  existence,  some  knowledge  about  the  inverse  operators  of 
M  e  2R  is  sufficient. 
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THEOREM  .3:  Let  the  following  assumptions  be  satisfied. 

1.  a)  All  operators  M  z  2W  have  an  inverse. 

b)  The  ranges  of  all  operators  M  z  2H  have  at  least  one 
element  r  ^0  in  common. 

c)  For  each  pair  MQ,  M^  of  operators  in  2 R  there  exists 

a  jamily  M(t)  e  2W  (0  <  t  <  1)  such  that  M(0)  =  M  , 

_  _  0 

M(l)  =  Mj  ana  such  that  the  family  M_1(t)r  depends 
continuously  on  t. 

At  least  one  operator  M  c  2)?  is  inverse-positive. 

Then,  all  operators  M  e  3)|  are  inverse-positive . 

This  theorem  is  an  immediate  consequence  of  the  lemma,  because 
the  elements  z(t)  =  M  (t)r  satisfy  all  required  assumptions. 

2.4  SPLIT  OPERATORS 

If  a  linear  operator  M,  occurring  in  some  given  problem,  is  not 
inverse-positive,  the  problem  can  often  be  transformed  such  that  it  is 
described  by  an  inverse-positive  operator  M. 

Define  R  to  be  the  linear  space  of  ordered  pairs  u  =  (u^u^  of 
elements  u^  u2  c  P>  with  the  order  relations: 
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u  ^  0 

iff 

u,  >  0, 

u_  <  0, 

1  — 

2  — 

u  y  0 

iff 

u!  y  °» 

u2*<  0; 

and  let  S  be  similarly  defined  using  ordered  pairs  of  elements 
of  S. 

Let  M^,  M2  denote  a  pair  of  linear  operators  mapping  P  into 
S  such  that 

M  =  -  M2,  (2.3) 

and  define  by 

Mu  =  (M^  -  M2u2,  M1u2  -  M2u1)  (2.4) 

an  operator  M  which  maps  R  into  S. 

Then,  one  may  ask  whether  this  operator  M  is  inverse-positive, 
that  means,  whether  for  u^,  u2 e  R, 

i  U1  -  ° 

imply  J  (2.5) 

I  u2-n' 

To  prove  this  property,  one  may  apply  Theorem  1.  The  following 

assumption  is  sufficient  for  Assumption  II  applied  to  M: 

Assumption  II’:  There  exists  an  element  z  e  P.  such  that 
z  1  0,  (Mx  +  M2)z  >»  0. 

This  follows  from  the  fact  that  z  =  (z,  -  z)  satisfies 
Assumption  II  applied  to  M. 
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SPECIAL  CASE  2:  Let 

M  -  A  -  B 

with  linear  operators  A,B,  and  suppose  that  A  and  B  can  be 
split  into 

A  ■  A^  -  A^t  B  ■  -  B^, 

such  that  A  defined  by 

Au  *  (Aiui  "  A2u2*  A1u2  “  A2U1^ 

satisfies  Assumption  I, and  the  operators  are  Posit^ve-  Then 

M  defined  by  (2.3)  with 

Ml-Al-V  M2  ~  A2  “  B2 

satisfies  Assumption  I. 

For  example,  if  A  itself  satisfies  Assumption  I,  then  A 
constructed  with  A^  ■  A,  A^  x  0  satisfies  Assumption  I. 

2.5  ERROR  ESTIMATION 

Let  be  given  an  equation 

Mu  =  r  (2.6) 

which  has  a  solution  u* e  R,  and  suppose  that  <f>  is  an  approximation 
of  u*  with  defect 

d[<f>]  *  -M<f  +  r. 

ERROR  ESTIMATION  1:  If  M  is  inverse-positive,  construct  an  element 
z  >  0  in  R  and  calculate  numbers  A,y  such  that  -A  <  y  and 
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-XMz  <_  d(t]  <.  uMz.  (2.7) 

Then ,  the  error  estimation 

-Xz  <_  U*  -  <J>  <_  uz 

holds. 

One  proves  this  statement  using  d[$]  -  M(u*  -  $) . 

REMARK:  The  element  z  has  to  satisfy  the  inequalities 
z  ^  0,  Mz  >  0 

which  are  very  similar  to  (2.1)  (except  in  the  trivial  case  d(<f]  =  0). 
Very  often,  the  element  z  used  in  the  error  estimation  will  satisfy 
(2.1).  In  other  cases,  an  element  z  satisfying  Assumption  II  can  be 
gained  by  slightly  modifying  the  element  z  used  in  the  error  estimation. 

If  M  is  not  inverse-positive,  one  may  get  an  error  estimation  using 
the  method  of  splitting  of  Section  2.4.  Let  M  be  an  operator  as 
described  in  that  section,  and  let  M  be  inverse-positive.  Then,  the 
only  solution  of 

Mu  =  r  with  r  =  (r,r) 

is  u*  =  (u*,u*). 

ERROR  ESTIMATION  2:  Let  (2.5)  be  satisfied  for  a  certain  splitting 
(2.3).  Construct  an  element  z  >_  0  in  R  and  calculate  numbers  X,u 
such  that  -X  <  u  and 


“X(M^  +  >l^)z  <_  d [ <i> ]  <_  p(M^  +  M2)z. 


(2.8) 
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Then t  the  following  error  estimation  holds , 

-Az  <_  u*  -  $  <_  uz. 

One  proves  this  statement  by  applying  Error  Estimation  1  to 

^  ^  A 

M,  4>  -  (♦,♦)»  Z  -  (z,-z). 


2.6  ITERATIVE  PROCEDURES,  EXISTENCE 

If  M  is  inverse-positive,  the  equation  Mu  ■  r  has  at  most 
one  solution.  Under  additional  assumptions,  one  can  prove  also  the 
existence  of  such  a  solution  using  the  means  of  the  preceding  sections. 

SPECIAL  CASE  1:  Consider  the  special  case  M  «  A  -  B  at  the  end 
of  Section  2.1  with  A  satisfying  Assumption  I  and  B  positive. 
Assume,  moreover,  that  M,  and  therefore  also  A,  satisfies 
Assumption  II  and  let  AR  ■  S.  Then,  A  is  inverse-positive,  and  the 
operator  T  defined  by  Tu  *  A  ^r  +  A  ^Bu  is  an  isotonic  mapping  of 
R  into  R. 

Suppose,  that  the  inequalities  (2.7)  in  Error  Estimation  1  are 
satisfied  and  define  sequences  {x^},  {y^}  by 


x(0)  ■  *  •  Xz-  *(0)  ■  *  +  « 

Ax(n+1)  ’  r  +  Bx(n)’  Ay(n+1)  '  r  +  By(n)  (n  '  0,1,2,...). 


i  .e . 


X(n+1)  *  Tx(n)’  y(n+l)  "  Ty(n) 


(n  -  0,1,2,...).  (2.9) 
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Then,  x^  ^y^,  and  by  applying  A  1  to  (2.7),  one  gets 
^x^^,  y(l)  —  y(0)‘  Because  T  is  isotonic,  these 
inequalities  yield 


X(0)  -  X(l)  - 


<  x.  .  < 

-  (n)  — 


i  y 


(n)  1 


i  y 


<  y /n\’  (2.10) 


(1)  -  y(0)* 


In  many  cases,  one  can  conclude  from  (2.9)  and  (2.10)  that  the 
sequences  {x^},  ^y(n)^  converge  to  fixed  points  x*  *  Tx*  and  y*  *  Ty*, 

respectively.  Each  fixed  point  is  a  solution 
of  the  given  equation  Mu  ■  r.  Because 

is  inverse-positive,  x*  =  y*.  (If  only  A,  instead  of  M, 
satisfies  Assumption  II,  we  cannot  conclude  x*  =  y*  which,  of  course, 
is  not  necessary  to  prove  the  existence  of  a  solution.) 

Under  suitable  additional  assumptions  on  the  space  R  and  the 
operators  A,B,  the  inequalities  (2.7)  of  Error  Estimation  1  yield 
the  existence  of  a  solution  u*  of  the  given  equation  Mu  =  r. 

SPECIAL  CASE  2:  Let  now  M  be  as  in  the  Special  Case  2  at  the  end 
of  Section  2.4.  Suppose,  moreover,  that  M  satisfies  Assumption  II' 
also,  and  let  AR  =  S.  (The  last  relation  holds,  for  example,  if  A  =  A^ 
and  AP  =  S.)  If  then  the  inequalities  (2.8)  of  Error  Estimation  2  are 
satisfied,  let  {x^},  (y^)  tBe  se9uences  defined  by 


x(0)  =  $  "  y(0)  =  <P  +  uz. 

AlX(n+l)  *  V(rH-l)  =  r  +  V(n)  *  B2y(n)’ 

Aly(n+1)  “  A2X(nfl)  =  r  +  Bly(n)  "  B2X(n) 
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These  sequences  again  satisfy  (2.10),  and,  therefore,  one  very  often 
can  conclude  that  these  sequences  converge  to  elements  x*,v*  with 

A^x*  -  A2V*  =  r  +  B^x*  -  B2y*. 

A^y*  -  A2X*  =  r  +  B^y*  -  B2X*. 

The  last  relations  are  equivalent  to 
Mx*  =  r  with  x*  =  (x*,y*),  as  well  as  to  My*  =  r  with 
y*  =  (y*,x*).  Because  M  is  inverse-positive,  x*  =  y*,  i.e.  x*  *  y*. 
Consequently,  u*  *  x*  =  y*  is  a  solution  of  the  given  equation.  (If 
only  A  =  A^  -  A2,  instead  of  M  =  -  M2,  satisfies  Assumption  II', 

one  cannot  conclude  x*  =  y*  but  still  can  prove  the  existence  of  a 
solution,  namely,  u*  =  *s(x*  +  y*).) 

Under  suitable  aaaitional  assumptions  on  the  space  R  ana  the 
operators  the  inequalities  (2.8)  of  Error  Estimation  2 

yield  the  existence  of  solution  u*  of  the  given  equation  Mu  =  r. 

3.  PROBLEMS  WITH  A  SI MOLE  DIFFERENTIAL  EOLATION  OF  AT  MOST  THE  SECOND  ORDER 
3.1  SUFFICIENT  CONDITIONS  FOR  INVERSE-POSITIVITY 

Consider  the  boundary  value  problem 

L[u]  =  L(u] (x)  =  a(x)u"  +  b(x)u’  +  c(x)u  =  s(x)  (0  <  x  <  1), 

V[uJ  =  au’(O)  +  Su(0)  =  A,  W[u]  =  yu'(l)  +  5u(l)  =  B. 

with  given  finite  real  functions  a,b,c,  and  s.  We  ask  for 
solutions  u  e  P  where  R  is  a  linear  subset  of  functions  in  C[0,1] 
which  have  sufficient  differentiability  properties.  At  least,  we  will 


(3.1) 
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require  that  u  e  R  possesses  all  derivatives  which  actually  occur 
in  the  given  problem.  (For  example,  u"(x)  shall  exist  if  a(x)  4  0 
for  that  point  x.) 


Under  what  conditions  do  the  inequalities 


L[u]  (x)  ■>  0  (0  <  x  <  1)  , 

V[u]  >  0,  W[u]  >  0 


imply  u(x)  ^  0  (0  x  <_  1)  (3.2) 


for  u  c  R? 


SPECIAL  CASE  A:  The  initial  value  problem 

L[u]  =  u'  +  c(x)u  *  s(x)  (0  <  x  <_  1) , 

V [u]  =  u(0)  =  A 

can  be  written  in  the  form  (3.1)  with 


W[u]  ■  L[u]  (1) ,  B  *  s(l) . 

SPECIAL  CASE  B:  The  singular  boundary  value  problem 
L[u](x)  =  a(x)u"  +  b(x)u*  +  c(x)u  *  s(x)  (0  <  x  <  1) 


with 


a(0)  =  a(l)  =  0, 


but  without  given  boundary  conditions,  can  be  written  in  the  form  (3.1) 
by  defining 

V [ u  J  =  L[u] (0) ,  A  =  s (0) ;  W[u]  -  L[u](l),  B  -  s(l). 
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THEOREM  4 :  If 

a(x)  <_  0  (0  <  x  <  1) , 

a  <_  0,  y  >_  0 

and  if  there  exists  a  function  z  e  R  such  that 

z(x)  1  0  (0<x<l), 

L[z](x)  >0  (0  1  x  <_  1),  V[z ]  >  0,  W[z]  >  0, 

then ,  (3.2)  holds  for  u  e  R. 


(3.4) 

(3.5) 


(3.6) 


PROOF:  Let  S  be  the  space  of  vectors  U  =  (s (x) ,  A,  B)  with 

real  functions  s(x) ,  defined  on  (0,1),  and  constants  A,B. 

Then,  the  operator  Mu  =  (L[u](x),  V[u],  W[u])  maps  R  into  S. 

Define,  moreover. 


u  >_  0  iff  u (x)  >_  0  (0  1  x  1)  ; 

u  >»0  iff  u(x)  >0  (0  <  x  <  1); 

0^0  iff  s  (x)  0  (0  <  x  <  1),  A  ^  0,  E  >_  0; 

U  ^»0  iff  s (x)  >0  (0  <  x  <  1),  A  >  0,  B  >  0. 


Then,  (3.2)  holds  iff  M  is  inverse-positive,  and  it  suffices  to 
prove  that  M  satisfies  the  Assumptions  I  and  II  of  Theorem  1.  To 
prove  Assumption  I,  let  us  suppose  that  w  ^  0,  Mw^-0,  but  not 

w  ^0.  Then,  w(x)  assumes  the  minimum  w(x^)  =0  at  some  point 
x0  *-  [°.l]  •  If  xqE  (0,1),  then  w'(x0)  =  0,  w"(x0>  1  0,  and 

therefore,  L[w](x0)  ^  0.  If  xQ  =  0,  then  w’(0)  >_  0  ,  and  therefore, 
V[w]  <_  0.  If  x^  =  1,  then  w'(l)  <_  0 ,  and  therefore,  W[w]  0. 
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Thus,  in  each  case  one  gets  a  contradiction  to  Mw  ^  0. 

For  the  SPECIAL  CASE  A,  the  assumptions  (3.4),  (3.5)  are 
always  satisfied.  Assumption  II  is  equivalent  to  (3.6). 

For  the  SPECIAL  CASE  B,  the  assumption  (3.5)  is  equivalent  to 

b(0)  <_  0,  b (1)  >_  0. 

For  example,  these  inequalities  hold  for 

L[u]  =  -x(l-x)u"  +  o(x-s5)u'  +  c(x)u, 

if  a  is  a  non-negative  constant.  If,  however,  the  constant  o  is 
negative,  one  needs  additional  boundary  conditions  in  order  to  form 
an  inverse-positive  operator  M. 

EXAMPLES  of  functions  z:  The  following  functions  z  satisfy 
the  inequalities  (3.6)  if  the  corresponding  conditions,  stated  below, 
hold . 

1)  z(x)  £  1, 
if 

c (x)  >0  (0  <  x  <  1),  B  >  0,  6  >  0;  (3.7) 

2)  z(x)  =  cos(Ti-e)  (x-b) 


with  small  enough  c  >  0,  if 
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2 

L[u]  =  -u"  +  c(x)u  with  c(x)  >  -tt  (0  <  x  <  1), 
V[u]  =  u(0),  W[u]  =  u(l); 


3) 


z(x) 


Nx 

e 


(3.8) 


in  case  of  the  initial  value  problem  (3.3)  with  c(x)  >  -N  (0  x  1). 

Further  functions  z  which  can  be  used  under  weaker  conditions 
as,  for  example,  (3.7)  are  given  in  [6]. 


3.2  CONNECTED  CLASSES  OF  EQUATIONS 

In  this  Section  3.2,  we  consider  only  operators  L  with  coefficients 
a(x),  b(x),  c(x)  which  are  defined  and  continuous  on  [0,1].  Moreover, 
we  restrict  ourselves  to  boundary  value  problems  (3.1)  with 

a(x)  =  -1  (0  <  1),  a  0,  y  >_  0,  (3.9) 

and  we  define  R  =  C^[0,1]  flC2(0,l). 

Each  triple  M  =  (L,  V,  W)  is  then  described  by  its  coefficient 
vector 


f  =  (b(x) ,  c  (x) ,  a,  6,  Y »  <5)- 

Let  9)1  be  a  set  of  suoh  triples  M  =  (L,V,W)  such  that  the 
corresponding  coefficient  vectors  f  form  a  connected  set  ,  i.e. 
for  each  pair  fQ,  e  f$,  there  exists  a  (component-wise)  continuous 
curve  f(t)  e  $(0  ^  t  1)  with  f  0  =  f  (0) ,  f^  =  f(l).  For 


20 


example,  9JI  may  be  a  convex  set. 

THEOREM  5:  Let  the  following  assumptions  be  satisfied. 

1)  For  each  M  =  (L,V,W)  r  2W*  the  homogeneous  problem, 
corresponding  to  (3.1),  has  no  nontrivial  solution. 

?,)  There  is  at  least  one  M  =  (L,V,W)  e  2)?  such  that  (3.8)  is 

satisfied. 

Then,  the  implication  (3.2)  is  true  for  all  M  *  (L,V,W)  e  2)?. 

PROOF:  The  statement  of  this  theorem  can  be  derived  from  Theorem  3 
by  considering  each  triple  M  =  (L,V,W)  as  an  operator  defined  as  in 
the  proof  of  Theorem  4,  and  using  the  following  facts. 

Because  of  (3.9)  each  M  e  9JI  satisfies  Assumption  I  according  to 
the  proof  of  Theorem  4.  Because  of  Assumption  1,  all  operators  M  c  2W 
have  an  inverse  M  ^  which  can  be  applied  to  r  *  (s(x),  A,B)  with 
s ( x )  "1,  A  =  1,  B  -  1. 

Moreover,  for  any  family  of  operators  M(t)  e  SD?  (0  <_  t  <_  1)  with 
continuous  coefficient  vectors  f(t),  the  functions  z(t)  =  M  \t)r 
depend  continuously  on  t. 

Finally,  at  least  one  M  r  2)?  is  inverse-positive,  because  (3.8) 
holds  for  at  least  one  K  r  ®l,  and  for  this  operator  M  Assumption  II 
is  satisfied  also. 

Roughly  spoken,  Theorem  5  says:  If  one  starts  with  a  triple 
(L,V,un  such  that  (3.2)  holds,  and  then  changes  the  coefficients 
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continuously,  the  implication  (3.2)  remains  true  as  long  as  one 
does  not  hit  a  triple  (L,V,W)  such  that  the  corresponding  problem 
(3.1)  does  not  have  a  unique  solution. 

A  very  simple  EXAMPLE  is  the  set  SDZ  of  triples  M  =  (L,V,W) 
with  (3.8). 

3.3  ERROR  ESTIMATION 

Consider  a  problem  (3.1)  such  that  (3.2)  holds.  Suppose  that 
there  exists  a  solution  u*  e  R,  and  let.  <J>  e  R  be  an  approximate 
solution.  Then,  one  can  get  a  bound  for  the  error  u*  -  $  in  the 
following  way. 

ERROR  ESTIMATION:  Construct  a  function  z  c  R  and  calculate 
a  constant  A  such  that 

|-L[<J)](x)  +  s(x)|  <_  AL[z](x)  (0  <  x  <  1), 

I  “V  [  <f> )  +  A  |  <_  A  V  [  z  ] 

|-W[<f>]  +  B  |  <  AW[z ]  . 

Then ,  the  error  estimation 

|u*(x)  -  (x)  |  <_  Az(x)  (0  1  x  <_  1)  (3.10) 

holds. 

3. A  A  DIFFERENT  APPROACH 

For  ordinary  differential  equations,  one  often  will  choose  an 
approximate  solution  <p  which  satisfies  the  given  boundary  conditions. 
Then,  one  does  not  need  (3.2)  but  only  the  weaker  property: 
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L[u](x)  ^  0  (0  i  x  i  1)  I 

I  imply  u(x)  _>  0  (0  <_  x  <_  1)  (3.11) 

V(u]  *  Wju]  ■  0  I 

for  a  e  R.  (We  now  suppose  that  the  coefficients  a(x),  b(x),  c(x), 
and  s(x)  are  defined  on  [0,1].) 

THEOREM  6:  If 

a(x)  <^0  (0  <_x  <_  1), 

if  thi  re  exists  a  function  z  e  R  such  that 

z(x)  ^  0  (0  f  x  ±  1) , 

L[z ]  (x)  >0  (0  1  x  <_  1),  V[z]  -  W[z]  -  0, 

(3. LI)  holds  for  u  e  R. 

This  theorem  can  be  proved  in  a  similar  way  as  Theorem  4.  One 
n  v  defines  Mu  *  L[u]  on  the  space  R  of  functions  u  e  R,  which 
sit isfy  the  homogeneous  boundary  conditions,  and  then  applies  Theorem  1 
(replacing  there  R  by  R) . 

Notice,  that  for  (3.11),  the  conditions  (3.5)  are  not  required. 

Using  this  different  approach,  one  can  consider  CONNECTED  CLASSES 
oi  equations,  similarly  as  in  Section  3.2,  and  one  can  also  derive  a 
corresponding  ERROR  ESTIMATION: 

| -L(4>] (x)  +  s(x)  |  <_  AL[z](x)  (0  1  x  <_  1) 

'••■plies  ju*(x)  -  (x)  |  <_  Az(x)  (0  f.  x  1)  , 

provided  $  e  R  satisfies  the  given  (inhomogeneous)  boundary  conditions 
and  z  c  R  satisfies  the  corresponding  homogeneous  boundary  conditions. 
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3.5  ESTIMATION  OF  THE  DERIVATIVE  OF  A  SOLUTION 

It  is  often  possible  to  calculate  bounds  for  the  derivative  of 
the  solution  of  a  problem  (3.1).  For  illustration,  we  consider  the 
special  case 

L[u](x)  =  -u"  +  b(x)u'  +  c(x)u  *  s(x)  (0  x  <  1) 
u(0)  *  A,  u(l)  =  B. 

(A  more  general  discussion  will  be  given  elsewhere .) 

Suppose  there  exists  a  solution  u*  e  R  *  C^[0,1]  fl  C2(0,l)  and 
there  has  been  gained  an  error  estimation 

|u*(x)  -  (J>(x)|  Az(x)  (0  x  <_  1)  (3.12) 

for  an  approximation  4  e  R  for  u*.  Let  z  e  R  satisfy  the 
homogeneous  boundary  conditions  z(0)  =  z(l)  =  0.  Let,  moreover, 
b(x)  be  bounded:  |b(x)|  <_  N  (0  f  x  1  !)•  Then,  one  can  calculate 
a  bound  for  (u*  -  $) ' . 

ESTIMATION  OF  THE  DERIVATIVE:  Construct  a  function  ;  e  0^0,1] 
and  calculate  a  constant  p  >_  0  such  that 

I -L [ cf) ]  (x)  +  s (x)  |  <_  nU'(x)  -  b(xK(x)]  -  A  | c (x)  | z (x)  (0<  x  _<  1),  (3.13) 

Az'(O)  <  p;(0)  .  (3.14) 

Then , 

|u*’  (x)  -  <j>'  (x)  I  ^  p;(x)  (0  ^  x  1  1) . 

PROOF:  Because  of  (3.12),  (3.13),  (3.14),  we  get  the  following 
estimations  for  v*  =  u*': 
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(v*  -  $')'  -  b(v*  -  p')  =  L  [  ^  ]  -  s  +  c(u*  -  p) 
u(r, '  -  hr)  (0  1  x  <_  1) , 

(v*  -  <J>’)(0)  *  lim  li  ^(u*(h)  -  tf> ( h ) ) 
h  ■*  +0 

<_  A  lim  li  ^z(h)  =  Az’(0)  <_  pr,(0) 
li  -»  +  0 

According  to  Theorem  4,  the  implication  (3.2)  holds  for  the 
operators 

L'[v]  =  v'  -  bv,  defined  on  (0,1), 

V ' [ v )  =  v(0),  W * [ v )  =  L'[vJ(l) 

(see  Special  Case  A  and  notice  that  Assumption  II  is  satisfied  with 
z  =  e‘  .)  Therefore,  the  inequalities  proved  above  yield 
v*  -  P*  <_  u4  (0  1  x  1).  In  a  similar  way,  one  shows  that 
-uf,  <_  v*  -  p  '  (0  x  <_  1)  which  then  proves  the  error  estimation. 

REMARK:  The  conditions  (3.13),  (3.14)  may  be  replaced  by  similar 
conditions  which  contain  an  inequality  at  x  =  1,  instead  of  (3.14). 

4.  SYSTEMS  OF  EQUATIONS  OF  AT  MOST  THE  SECOND  ORDER 
4.1  SUFFICIENT  CONDITIONS  FOR  INVFPSE-POSITIVITY 

We  consider  now  a  boundary  value  problem  for  m  unknown  functions 

u  (x) , . . . ,u  (x)  : 
i  m 
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LjM  =  a^x)^’  +  b.(x)uj  +2  cik(x)uk  =  Si(x)  (0  <  x  <  D» 

k-1 

m 

V"1  -  V;(0>  +  V  6ikuk(0)  -  Alt 

k*  i 

m 

Wilu]  =  Vi(1)  +  2  6ikUk(1)  =  Bi  <i  -  1.2 . m). 

k-1 

Let  R  be  a  linear  set  of  m-tuples  u  -  (u  (x) ,  u,(x) . u  (x)) 

*■  *>  m 

such  that  all  uk  e  C[0,1]  and  such  that  for  u  e  R,  the  components 
uk(x)  possess  all  derivatives  which  actually  occur  in  the  given 
problem.  We  ask  for  solutions  u  e  P. 

Again,  this  problem  includes  initial  value  problems  and  singular 
boundary  value  problems  as  special  cases. 

THEOREM  4':  Suppose  that 

ai<x>  1°  (0  <  x  <  1),  ^  <  0,  yj  >  0  (4.1) 

end 

cik*x'  1  0  <  x  <  D.  8lk  £  0,  6lk  ^  0  for  i  4  k  (4.2) 

(i.k  =  1,2, .. . ,m) . 

Moreover,  let  there  exist  z  -  (z ...... z  )  e  R  such  that 

i  m 

zf(x)  ^  0  (0  1  x  £  1), 

L1[z](x)>0  (0  <  x  <  1),  V A [ z J  >  0 ,  Wt[zJ  >  0,  (i  =  1,2,. ..,»). 
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Then ,  the  inequalities 

L^[u](x)  ^  0  (0  x  1) 

v± Iu ]  1  0,  l^fu]  _>  0  (i  -  1,2,. ...m) 

imply 

u^(x)  >_  0  (0  £  x  <  1;  i  *  l,2,...,m). 

This  theorem  can  be  proved  by  generalizing  the  proof  of  Theorem  4. 
One  defines  M  to  be  the  operator 

Mu  =  (L. [u] , . . . ,L  [u] ;  V  [u] , . . . ,V  [u] ;  W(u] , . . . ,W( u] ) 
i  ml  ml  m 

(A.  3) 

*  (Li(u],  V1(u],  W.[u]), 

which  maps  R  into  the  space  S  of  elements  U  *  (s^(x),  A^,  B^) 

(i  *  l,2,...,m)  with  functions  Sj(x)  defined  on  (0,1)  and 
constants  A^,  B^.  The  inequalities  u  ^  0,  U  ^  0  shall  be  understood 
as  holding  component-and  pointwise. 


4.2  SPLITTING  OF  THE  PROBLEM 

If  the  assumptions  of  Theorem  4'  are  satisfied,  one  can  get  an  error 
estimation  for  an  approximate  solution  $.  In  many  cases,  however,  an 
error  estimation  is  also  possible,  if  (4.2)  does  not  hold.  We  apply  the 
method  of  splitting,  described  in  Section  2.4,  to  the  operator  M  in 
(4.3).  Let 
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Wu[u]  -  v.ujd)  ♦  6ilui(l)  -  I ’  5jk  uk(l). 


Li2  “  Lil  “  Lf 


12 


11 


V  w  «  w  -  w 

V.  wi2  Wil 


(1  m  1,2,... ,m) 


where  we  have  used  the  notations 
m 

I’  *2’  f”  *  *2(lfl  -  f>  * 

k=l 

kyi 

Define  then,  similar  as  in  (4.3), 

M)U  =  (Ln  tu) ,  Viilul»  W.1(u]  ),  M2u  =  (L.2[u],  Vi2[u],  k’i2  [u] ) . 

Then,  for  the  operator  M,  corresponding  to  the  splitting  M  *  -  M2 

(see  (2.4)),  Assumption  I  is  always  satisfied  if  only  (4.1)  holds. 

The  corresponding  operator  M  =  +  M2  is 

Mu  *  (Lilu],  V^[u] ,  VJ  ^  l  u  ] ) 

with 


L.[u]  =  aiU”  +  t^u’  +  cuu.  -  S’|cik|uk, 

V^u]  =  0^(0)  +  Bj.u.W  -  I'  |  Bjkluk(°)  ♦ 

Wjtu]  -  Y.^d)  +  6^(1)  -  Z’|6lk|uk(l). 

With  these  notations.  Assumption  II'  is  as  follows: 

Assumption  II':  Let  there  exist  z  =  (z^,...,zm)  e  P  such  that 

Zi(x)  >  0  (0  1  x  <  1), 

L. [z] (x)  >0  (0  <  x  <  1),  V.[z]  >  0,  W± [ z ]  >0  (i  =  1,2,... ,m) . 
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4.3  EPROR  ESTIMATION' 

Let  the  given  proolem  in  Section  4.1  have  a  solution  u*  e  R, 
and  let  4>  e  R  denote  an  approximation  for  u*.  Suppose,  that  the 
inequalities  (4.1)  hold  and  that  Assumption  II'  in  Section  4.2  is 
satisfied.  Then,  one  can  get  a  bound  for  the  error  u*  -  <J>  in  the 
following  way. 

ERROR  ESTIMATION:  Construct  an  element  z  e  R  and  calculate  a 

constant  A  such  that 

|-Li[4»](x)  +  s.(x)|  <_  XLi(z)(x)  (0  <  x  <  1), 
l-V.m  +  Aii  1  AV^z], 

|-W±m  +  Bi|  <  AW^z]  (i  =  1,2,...  ,m)  . 

Then , 

| u  (x)  -  <j)i(x)  |  <_  Azi(x)  (0  <  x  <  1;  i  *  1,2, . . .  ,m) . 

For  z  in  Assumption  II',  such  a  constant  X  always  exists. 

This  statement  is  a  consequence  of  Error  Estimation  2  in  Section  2.5. 

5  BOUNDARY  VALUE  PROBLEMS  OF  THE  FOURTH  ORDER 
5.1  THE  GENERAL  PROBLEM 

Consider  a  differential  equation 

L [ u ] (x)  =  s(x)  (0  <  x  <  1) , 


together  with  two  boundary  conditions  at  x  =  0: 
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Vjlu]  =  A.  (i  =  1,2), 
and  two  such  conditions  at  x  =  1: 

VJ±[u]  -  Bi  (i  =  1,2). 

Let  L  be  a  linear  differential  operator  of  the  fourth  order  such 
that  the  coefficient  of  u ^  is  positive  in  [0,1],  and  let 
V  [u],  W^[u]  be  linear  combinations  of  derivatives  of  u  up  to  the 
third  order.  For  simplicity,  we  assume  that  the  coefficient-functions 
of  L  are  analytic  on  [0,1]. 

Under  uhat  conditions  do  the  relations 

L[u](x)  ^  0  (0  x  1) 

imply  u(x)  ^  0  (0  ^  x  <_  1)  (5.1) 

Vt[u]  -  Wt[u]  -  0  (i  *  1,2) 

for  all  analytic  functions  u  on  [0,1]  .' 

This  is  true  iff  the  corresponding  Green’s  function  exists  and  is 
non-negative.  Therefore,  if  this  implication  is  true  for  all  analytic 
functions  u,  it  is  also  true  for  all  u  e  [0,1] . 

There  are  some  difficulties  in  applying  Theorem  1  to  this  problem 
directly.  In  applying  Theorem  1  to  second  order  problems,  we  used  the 
fact  that  the  first  and  second  derivatives  of  a  function  have  to 
satisfy  certain  necessary  conditions  at  a  point  where  this  function 
assumes  an  extremum.  However,  at  such  a  point  nothing  can  be  said  about 
higher  derivatives  which  occur  in  problems  of  higher  order.  Therefore, 
we  split  the  given  operator  L. 
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We  construct  linear  differential  operators  of  the  second  order 
L^,L0  such  that  for  all  analytic  functions  u: 

L2  [v]  -  w(x)(L[u]  +  q(x)u)  (°  1  x  <_  1)  (5.2) 

with 


L1 [u]  -  v 


(5.3) 


and  some  fixed  (analytic)  functions  w,q. 


Then,  we  ask  for  conditions  such  that  for  all  analytic  u(x) : 


L[u] (x)  >0  (0  <  x  <  1) , 


V± [ u ]  =  W1(u]  *  0  (i  *  1,2) 


imply 


u(x)  >  0, 


v(x)  >  0  (0  <  x  <  1). 


(5.4) 


In  order  to  gain  such  conditions,  we  formulate  this  problem  in 
abstract  terms  and  then  apply  Theorem  1. 

Let  R  be  the  space  of  all  functions  u  which  are  analytic  on 
[0,1]  and  satisfy  the  boundary  conditions  V^[u]  =  W^[u]  *  0  (i  =  1,2), 
and  let  S  be  the  set  of  analytic  functions  U  on  [0,1],  Define 
for  u  e  R,  respectively  U  c  S: 


u  >  0  iff 


with  (5.3), 


u(x)  >_  0, 

v(x)  >0  (0  <  x  <  1) 


(5.5) 


U  10  iff  U(x)  _>  0  (0  <.  x  <  1) , 


and  let  u^O,  U  denote  the  corresponding  strong  order  relations, 

defined  by  Property  o  .  Then,  Mu  =  L[u]  is  a  mapping  of  R  into  S, 
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and  (5.4)  holds  Iff  M  is  inverse-positive. 

Thus,  we  may  apply  Theorem  1  in  order  to  get  sufficient  conditions 
for  the  property  (5.4).  It  is  not  possible  to  present  here  a  general 
theory  of  this  kind.  A  first  paper  on  this  topic  will  soon  be  published 
elsewhere  [9].  Instead  of  reviewing  general  results,  we  will  discuss 
here  a  special,  but  typical  case. 

5.2  A  SPECIAL  PROBLEM 

For  the  special  boundary  value  problem 

L[u](x)  =  uIV  -  bu"  +  c(x)u  *  s(x)  (0  x  <_  1) ,  (5.6) 

u(0)  -  A^  u'(0)  -  A2,  u(1)  -  B1,  u'(1)  -  B2  (5.7) 

with  constant  b,  let  us  choose  the  following  splitting: 

L. [u]  =  -pu"  +  p'u*  +  (P  +  y)u, 

(5.8) 

L2[v]  =  -pv"  +  p'v’  +  (P  -  y)v 

where  p(x)  denotes  an  analytic  function,  P  *  -bp"  +  bbp,  and  y 

is  a  constant.  For  these  operators,  (5.2)  (5.3)  hold  with 

2  2 
w  -  p  ,  q  =  L2 [P  +  y]  -  cp  . 

In  particular,  we  choose 

x 

p(x)  =  fgiOdf,  (5.9) 

0 

where  g  denotes  a  function,  such  that  for  some  constant  c^: 

TV  2 

8  -  2bg"  +  b  g  =  4cQg  (0  <_  x  <_  1)  (5.10) 


32 


g(0)  »  0,  g'(0)  *  1,  gCx)--g(l-x)  (0  <  x  <  1), 

g(x)  >0  (0  <  x  <  *s). 

Moreover,  let 

Y  -  -  *s(l-e)  (5.12) 

with  some  fixed  e  such  that  0  <  e  <  2. 


(5.11) 


THEOPEM  6:  Suppose  the  following  assumptions  are  satisfied. 

I.  For  some  constant  cQ  with 

c(x)  £  c0  !°  i  *  i  1).  (5.13) 

there  exists  a  solution  g  of  the  equation  (5.10)  satisfying  (5.11). 
Let  then  L^,L9  denote  the  operators  described  above. 

II.  There  exists  a  function  z  e  R  such  that 


L[z]  (x)  >  0  (0  £  x  <_  1) ,  I 

z(x)  >_  0,  [z]  (x)  >_  0  (0  <  x  <  1).  I 


(5.14) 


Then,  for  all  analytic  functions  u  on  [0,1],  the  relations 
IV 

L[u]  =  u  -  bu”  +  c(x)u  ^  0  (0  ±  x  f  1), 

u(0)  =  u'(0)  *  u(l)  =  u’(l)  =  0 

|u(x)  >  0,  (5.15) 

L-Ju]  (x)  >_  0  (0  1  x  <_  1) 


PROOF:  We  will  show  that  the  operator  M  defined  in  Section  5.1 
satisfies  the  assumptions  of  Theorem  1,  for  the  special  case  which  we 


consider  here. 
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The  functions  v  *  L^[u]  with  u  e  R  satisfy  v(O)  *  v‘(0)  = 
v ( 1 )  *  v'(l)  ■  0.  There  are,  however,  functions  v  of  that  kind  such 
that  v"(0)  t  0,  v"(l)  t  0.  Therefore,  and  because  u  e  R  satisfies 
the  homogeneous  boundary  conditions,  corresponding  to  (5.7),  for 
u  c  R: 

u(x)  >  0,  v(x)  >0  (0  <  x  <  1), 

u^-0  iff  (5.16) 

u"(0)  >  0,  u"(l)  >  0,  v"(0)  >  0,  vM(l)  >  0. 

Moreover,  for  U  e  R: 

U  ^0  iff  U(x)  >0  (0  <_  x  <_  1) . 

Thus,  Assumption  II  of  Theorem  1  is  satisfied  because  of  the 
conditions  (5.14),  and  Assumption  I  of  Theorem  1  is  equivalent  to  the 
following  condition:  For  u  e  R,  the  inequalities 

L[u]  (x)  >  0  (0  <  x  <_  1) 

(5.17) 

u (x)  21  °i  v(x)  =  L^[u](x)  0  (0  1  x  1) 

imply  the  inequalities  in  (5.16). 

The  function  p  in  (5.9)  is  positive  in  the  open  interval  (0,1) 

2 

because  of  the  properties'  of  g.  The  function  q^  =  L^fP  +  y]  -  c^p 
is  constant  because 

TV  ? 

2qJ}  =  p[g  -  2bg"  +  (b  -  4cQ)g]  5  0  (0  <  x  <  1). 

The  constant  value  of  qQ  is 

qo(x)  =  qo(0)  =  k[l  "  (1_c)2]  "  °‘ 
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Therefore,  and  because  of  (3.13), 

q  (x)  =  q0(x)  +  (c0  -  c(x))p2(x)  0  (0  I  x  <_  1). 

Let  now  the  inequalities  (5.17)  be  satisfied  for  some  u  c  R. 

Then,  the  function  on  the  right-hand  side  of  the  equality  (5.2)  is 
positive  in  (0,1),  and  therefore,  one  can  prove,  by  the  same  method 
that  was  used  for  second  order  equations,  that  v(x)  >0  (O  <  x  <  1) . 
But  then,  the  same  proof  applied  to  equation  (5.3)  yields 
u(x)  >0  (0  <  x  <  1). 

In  order  to  prove  u"(0)  >  0  and  v"(0)  >  0,  let  us  assume  that 
v"(0)  >0  is  not  satisfied.  Then,  v"(0)  =  0  where 

v"  (0)  =  (+  V'(0)  +  y)  u"  (0)  =  ru"(0),  (5.18) 

and  therefore,  u"(0)  =  0.  These  equalities  then  yield 

0  1  v'"  (0)  =  (-— p"(0)  +  y ) u"'  (0)  =  (-1  +  7)u"'(0)  and  u"'(0)  >0 

which  is  only  possible  if  v"'(0)  =  u'"(0)  =  0.  In  this  case,  however, 

TV  S  TV  r 

0  ^  v  (0)  =  (-jp"(0)  +  y ) u  (0)  =  (-3  +  y)L[u](0) 

which  contradicts  (5.17).  Therefore,  v”(0)  >  0  and  also  u"(0)  >  0 
because  of  (5.18). 

In  the  same  wav,  one  proves  u"(l)  >  0,  v"(l)  >  0. 

EXAMPLE  1:  For 

4c0  =  (4~2  +b)2, 

the  problem  (5.10),  (5.11)  has  the  solution 

R ( x )  =  -  ^7  sin  2 n  (x-?5)  . 
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EXAMPLE  2:  For  cQ  «  b  -  0,  the  problem  (5.10),  (5.11)  has 
the  solution  g(x)  ■  x(l-x) (l-2x) ,  so  that 

p(x)  -  %x2(l-x)2. 

2  2 

Moreover,  the  function  z  ■  x  (1-x)  satisfies  (5.14)  if  c(x)  =  0. 

5.3  CONNECTED  CLASSES  OF  PROBLEMS 

Theorem  3  cannot  be  applied  immediately  because  the  order  relations 
defined  in  R  now  depend  on  the  operator  M.  However,  one  still  can  apply 
the  basic  idea  of  that  theorem. 

Again,  we  consider  the  problem  (5.6),  (5.7);  but  first  we  assume 
that  c(x)  =  cq  *  const. 

Let  St  be  a  connected  set  of  coefficient  vectors  f  *  (b.c^), 
and  let  ?  denote  the  set  of  operators  L  with  (b,cQ)  e  . 

If,  in  the  following,  connected  with  some  operator  L  c  ?  ,  an 
operator  L^  is  considered,  this  operator  L^  shall  be  defined  by  (5.8) 
through  (5.12)  with  the  coefficients  (b,c^)  of  L.  For  example,  in 
(5.15)  this  operator  L^  is  to  be  used. 

THEOREM  7a:  Suppose  that  the  following  assumptions  are  satisfied. 

la)  For  each  L  c  f,  the  homogeneous  problem,  corresponding  to 
(5.6),  (5.7),  has  no  nontrivial  solution. 

lb)  For  each  L  e  2,  the  corresponding  problem  (5.10)  has  a 
unique  solution,  such  that  (5.11)  holds. 
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2.  For  at  least  one  L  e  9,  the  property  (5.15)  holds. 

Then,  the  property  (5.15)  holds  for  all  L  e  9. 

PROOF:  Suppose  (5.15)  does  not  hold  for  z  9  with  coefficient 

vector  I1,  but  does  hold  for  L°  e  9  with  1°.  Let  then 
ft  e  il  (0  <  t  <  1)  be  a  continuous  curve  which  connects  and 

f^,  and  let  Pt(x),  be  the  corresponding  quantities  defined  by 
(5.9)  through  (5.12).  Moreover,  define  z  R  to  be  the  solution  of 
L[z](x)  =1  (0  x  <_  1)  where  L  has  the  coefficient  vector  f^. 

The  functions  pt(x),  zt(x),  and  y  depend  continuously  on  t. 

Because  of  Theorem  6,  it  is  sufficient  to  show  that  the  functions 
z£(x)  satisfy  >  0  for  0<  x  <  1  in  the  sense  of  (5.5).  We  will 
even  show  that  z£  ^0  for  0  <_  t  <_  1 .  Notice  that  the  inequalities  on 
the  right-hand  side  of  (5.5),  respectively  (5.16),  partly  depend  on  L^. 
Of  course,  in  the  definition  of  zf  ^  0,  respectively  zt^0,  the 
corresponding  operator  constructed  with  f(t)  must  be  used. 

We  have  zQ  >_  0  because  (5.15)  holds  for  L  *  L°.  If  z^^O 

(0  i  t  i  1)  is  not  true,  then,  because  of  continuity  conditions,  there 

is  a  smallest  tfl  z  [0,1]  such  that  z  ^  0,  but  not  z  ^  0.  Such 

C0  C0 
a  number  tQ,  however,  cannot  exist  because  the  operator  M, 

corresponding  to  t0,  satisfies  Assumption  I  of  Theorem  1,  so  that 

z  ^  0  implies  z  X  0. 

C0  C0 

We  can  extend  the  result  of  Theorem  7a  to  operators  with  variable 
coefficient  c(x).  Suppose  that  the  assumptions  of  Theorem  7a  are 


satisfied. 
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THEOREM  7h:  Let  b  and  denote  fixed  numbers  such  that 

(b,cQ)  e  St.  If  the  homogeneous  problem,  corresponding  to  (5.6), 
(5.7),  has  no  nontrivial  solution  for  all  (b,c(x))  with 

C(x)  <  c (x)  c0  (0  <  x  <  1),  (5.19) 

where  C(x)  denotes  some  fixed  function,  then  for  all  operators  L 
with  such  coefficients  the  property  (5.15)  holds. 

The  proof  Is  similar  to  the  proof  of  Theorem  5.  Notice,  that 
now  Lj  is  a  fixed  operator. 

We  will  now  give  a  more  explicit  example.  Let,  for 

2 

-10w  <_  b  <  °°,  X  ■  A*(b)  denote  the  smallest  eigenvalue  of  the 

eigenvalue  problem 

TV  2 

g  -  2bg"  +  b  g  -  4Ag,  g(x)  *  -g(l-x)  (0  <  x  <  1) , 

(5.20) 

g(0)  -  g’(0)  -  g(l)  -  g’(l)  -  0, 

and  suppose  that  X  ■  Ai(b)  is  the  largest  eigenvalue  of  the 
problem 

uIV  -  bu"  *  -Xu  (0  <  x  <  1), 

(5.21) 

u(0)  -  u’(0)  -  u(l)  -  u’(l)  =  0. 

Then,  the  following  Corollary  is  a  consequence  of  Theorem  7 


(see  Figure  2). 
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COROLLARY:  If 
2 

-IOtt  <  b  <  00  ana  A^  (b)  <  c(x)  <_  A*(b)  (0<  *  <  1),  (5.22) 

then,  for  all  u  e  C^[0,1], 

L[u ] (x)  >  0  (0  ^  x  <  1), 

u(0)  =  u'(0)  *  a(l)  ■  u' (1)  «  0 

PROOF:  Let  St  denote  the  connected  set  of  vectors  f  *  (b,Cg) 

with 

2 

-10n  <  b  <  »  ,  Aj(b)  <  cQ,  0  <_  cQ  <  A*(b), 

and  let  8  be  the  set  of  the  corresponding  operators  L.  We  will  show 
that  the  assumptions  of  Theorem  7a  are  satisfied  for  this  set  8. 

Because  for  each  L  e  8*  A^(b)  <  c^,  the  homogeneous 

problem  corresponding  to  (5.6),  (5.7)  has  no  nontrivial  solution.  Thus, 
Assumption  la  of  Theorem  7a  is  satisfied. 

Because  c^  <  A*(b)  for  each  L  e  8  ,  the  corresponding  problem 

IV  2 

g  -  2bg"  +  b"  g  =  4cQg  (0  1  x  <_  1), 
g(0)  -  g(l)  =  0,  g’(0)  -  g’(l)  «  1 

has  a  unique  solution  g  which,  obviously,  satisfies  g(x)  *  -g(l-x). 
Therefore,  in  order  to  prove  that  Assumption  lb  is  satisfied,  we  have 
only  to  show  that 

g (x)  >0  (0  <  x  <  -2) .  (5.24) 


imply  u(x)  >  0  (0  <.  x  <_  1). 

f5.231 


This  inequality  is  true  for  the  function 


x(l-x)(l-2x)  which 
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FIGURE  2: 

Shaded  domain:  described  by  (5.25); 

-  :  eigencurve  A  =  cQ  «  A*(b)  of  problem  (5.20); 

-  .  eigencurves  A  =  cQ  «  A(b)  of  problem  (5.21). 
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corresponds  to  L  [u]  ■  u  .  Suppose,  the  inequality  does  not  hold 
for  some  g\  corresponding  to  an  operator  c  £  .  Then,  by 
continuity  arguments  of  a  similar  kind  as  we  have  used  before, 
one  can  conclude  that  there  exists  a  function  g  associated  with  an 
operator  L  e  ?  such  that  g(0  =  g'(0  *  0  for  some  £  e  (O,1^]. 

Suppose  first,  that  this  operator  L  has  coefficients  b,cQ 
2 

satisfying  b  <  4cQ.  Then, 

g(x)  ~  ~~qT  [sin  6  •  sinh  2a(x-^)  -  sinh  a  •  sin  2B(x-^)], 

where  2a  -  /  2/^  +  b,  26=  f  2/~c^  -  b,  and 

A(a, 8)  =  2a  sin  B*cosh  a  -  26  sinh  a*cos  6. 

The  value  A(a,B)  does  not  vanish  for  any  numbers  a,B  which 
belong  to  a  vector  (b,cQ)  e  ,  because  A(a,B)  =  0  determines  the 
eigenvalues  of  the  problem  (5.20),  and  cQ  <  A*(b). 

The  argument  £ ,  mentioned  above,  cannot  be  £  -  h  because 
g'(^)  =  0  would  yield  sinh  a/a  =  sin  B/B  which  is  not  satisfied 
for  any  a  >  0,  6  >  0. 

For  0  <  £  <  ht  we  calculate 

0  =  A(a,B)  [g(0  -Ba-cosh  2a(£-1i)  -  g'(0  •  sinh  2a(£-^)J 

=  -(£,-%)  sinh  a*A(a,B) 

with  a  ■  ta,  3  =  tB,  and  0  <  t  =  2(£-1s)  <  1.  However,  A(a,B)  f  0 
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because  the  vector  (b,cQ)  corresponding  to  a,B,  also  belongs 
to  JJ  . 

2 

This  proves  (5.24)  for  all  L  e?  with  b  <  4cQ.  The  other 
2 

cases,  where  b  >_  4Cq,  can  be  treated  similarly. 

Assumption  2  of  Theorem  7a  is  also  satisfied,  because,  as  we 

have  seen  at  the  end  of  Section  5.2,  the  property  (5.15)  holds  for 

IV 

the  operator  L[u]  =  u  which  belongs  to  ?. 

In  order  to  show  that  (5.15)  holds  for  all  u  e  R  if 
2 

-IOtt  <  b  <  °°  and  A^b)  <  c(x)  £  A*(b)  -  6(b)  (0  <_  x  <_  1)  (5.25) 


with  some  sufficiently  small  6(b)  >  0,  we  apply  Theorem  7b  with 

C(x)  =  Aj(b),  and  c^  =  A*(b)  -  6(b)  for  each  fixed 
2 

b  e  (-10ir  ,  °°) .  It  has  only  to  be  proved  that  the  equation 
XV 

u  “  bu"  +  c(x)u  =  0  yields  u(x)  =0  for  u  e  R.  This  can  be  done 
by  a  standard  procedure.  If  u(x)  f  0,  then 
1  1 


0  f  u(uIV  -  bu"  +  c(x)u)dx  =  f  (u"u"  +  bu'u'  +  c(x)u2)dx 
0  0 

>f  (u"u"  +  bu'u'  )dx  +  A^b)  J  u2(x)dx. 


0 


0 


This  is  a  contradiction  to  the  fact  that  A^b)  is  the  largest 
eigenvalue  of  the  problem  (5.21). 


Thus,  under  the  assumption  (5.25),  the  implications  (5.15)  and 
(^•23)  hold  for  all  u  c  R.  In  order  to  show  that  (5.23)  is  even  true 
for  all  u  e  C^[0,1]  if  the  weaker  assumption  (5.22)  is  satisfied,  one 
uses  the  fact  that  the  Green's  function,  corresponding  to  the  problem 
(5.6),  (5.7)  exists  for  c(x)  > 


A^(b)  and  depends  continuously  on 
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c(x),  and  that  this  function  is  non-negative  if  and  only  if 
(5.23)  is  true  for  all  u  e  R,  or  all  u  e  C^[0,1]. 
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6.  NUMERICAL  APPLICATION 

6.1  A  METHOD  OF  APPROXIMATION  AND  ESTIMATION 

Let  there  be  given  a  differential  equation 
L[u](x)  *  s(x)  (0<x<l) 

of  the  N-th  order  (for  example,  N  *  2,  or  N  *  4) ,  together  with 
appropriate  linear  (inhomogeneous)  boundary  conditions,  and  suppose 

that 

L[u](x)  ^  0  (0  x  1  1)  implies  u(x)  ^  0  (0  x  1) 

for  all  u  e  C^tO.l]  which  satisfy  the  corresponding  homogeneous 
boundary  conditions.  Assume,  moreover,  that  the  given  problem  has  a 
solution  u*  e  C^tO.l]. 

Then,  one  may  calculate  an  approximation  $  for  u*  and  bounds 
for  the  error  u*  -  <j>  in  the  following  way. 

STEP  A  (Approximation).  A  development 

4>  -  4>n  +  a  +  . . .  +  a  <p 
u  l  i  mm 

is  set  up  such  that  for  arbitrary  constants  a,,..., a  ,  the  function 

1  m 

<J>  is  contained  in  CN[0,1]  and  satisfies  the  given  boundary  conditions. 
The  constants  are  determined  such  that  the  defeat 

d[$]  (x)  =  -L[ <j>]  (x)  +  s(x) 

is  orthogonal  to  m  appropriately  chosen  functions  ^  (k  =  l,2,...,m), 
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with  respect  to  an  inner  product 


(u,v) 


/ 


w(x)u(x)v(x)dx. 


(6.1) 


where  w  e  C0[0,1]  denotes  some  weight  function.  That  means,  the 
constants  constitute  a  solution  of  the  linear  system 


a . 
1 


+ 


0 


(k  =  1,2, ... ,m) 


wi  th 


(6.2) 


8ik  =  =  -Lt  •i'i  J  (i»k  =  1*2  ... .  ,m)  , 

rk  =  ^“o’^k^*  “0  =  d^0]  = 

STEP  E  (Estimation).  A  function  z  e  C_.  [0,1]  satisfying  the 

ft 

homogeneous  boundary  conditions  is  chosen,  and  a  constant  X  >_  0 
is  calculated  such  that 


| cl [ <j) ]  (x)  |  <_  XL[z]  (x)  (0  x  <_  1)  .  (6.3) 

Then,  the  error  estimation 

|u*(x)  -  $ (x)  |  £  Xz (x)  (0  1  x  <_  1)  (6.4) 


holds . 

This  procedure  has  been  programmed  for  a  computer.  In  the  program, 
the  occurring  integrals  (6.1)  are  calculated  by  some  approximation 
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method,  and  the  system  (6.2)  is  solved  by  an  elimination  procedure. 

In  order  to  obtain  a  number  X  satisfying  (6.3),  the  function 
d[$]  (x)/L[z]  (x)  is  calculated,  and  the  corresponding  graph  is 
dotted  for  a  large  number  of  points  x. 

6.2  AN  EXAMPLE  OF  THE  SECOND  ORDER 
The  boundary  value  problem 

-u"  +  2xu'  +  (l-x2)u  =  1-x2  (0  <_x  <  1)  , 

u'(0)  -  u(l)  -  0 

has  the  exact  solution 

u*  *  1  -  exp^(x2-l)  . 

We  calculated  approximate  solutions  with 

i+1 

4>0  =  0»  i  *  1-x  (i  *  1,2,...  ,m)  , 

but  chose  w  and  the  in  three  different  ways  (Cases  1,2,3). 

For  the  error  estimation,  the  function  z  =  z°  with 

2  2 

z°(x)  *  ^(e  -  eX  )  and  L[z°]  -  eX  +  (l-x2)z°  (6.5) 

was  used.  The  Tables  1,2,3  contain  some  of  the  corresponding  results, 
namely 

6  =  |  ju*  -  $||,  and  X||z||  ■  4*e  with  X  -  | |d($]/L[z] | | , 
where  the  notation 

|  |  u  |  |  «  max  { |  u  (x)  |  :  0  <_  x  <_  1 } 

is  used.  Obviously,  x||z||  is  the  maximum  of  the  bound  for  the  error 
|u*  -  <p  | .  The  numbers  6  are  rounded  in  the  usual  way,  the  numbers 
X  J  | z  |  |  are  rounded  up  to  the  next  larger  decimal. 


The  simplest  weight  function  w(x)  =  1,  and  the  functions 
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k-1 

^  ■  x  (k  ■  l,2,...,m)  give  the  results  in  Table  1,  for 
m  ■  2, 3, 4, 5  (Case  1).  For  this  case,  the  maximal  values  of  the 
quotient  | d [  <f>]  (x)  |  /L[z ]  (x)  occur  near  the  boundary  points  x  *  0 
and  x  *  1. 

Therefore,  in  the  second  case,  the  same  functions  ij>^,  but 
weight  functions 

w(x)  -  [1  -  y(2x-l) 2 ]  1  (y  a  |-) 

are  used  which  have  larger  values  near  the  boundary  points  than 
in  the  middle  of  the  interval.  -The  corresponding  results  for  m  ■  4 
are  contained  in  Table  2. 

A  similar  effect  can  be  obtained  with  the  simple  weight  function 
w(x)  =  1,  but  rational  functions  t^.  For  m  =  4, 


<1*1  =  1»  =  x» 


and  different  values  of  c,  one  gets  the  results  in  Table  3. 

The  result  in  the  second  case,  for  y  ■  7/8,  comes  close  to  that 
which  one  would  gain  by  requiring  that  the  Tschebyschef f  norm  of 
d[<f>]  (x)/L[z]  (x)  is  minimal,  because  the  extrema  of  this  quotient 
have  alternate  signs  and  almost  equal  absolute  values. 

In  all  cases  which  we  treated,  we  obviously  did  not  get  very  much 
smaller  error  bounds  by  using  another  approach  than  the  simplest  one 


in  Case  1. 
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TABLE  1 

(w(x)  =  1,  i|>k  -  xk  1) 

MMI 

■  X • ■ 

m 

6  -  llu*  -  4>  1  1 

i  error 

bound  1  I 

3 

0.000  715 

0.033 

4 

0.000  069  0 

0.008 

3 

5 

0.000  005  04 

0.000 

92 

6 

0.000  000  525 

0.000 

15 

TABLE  2  (w  -  [1  -  y(2x  -  l)2]' 

-1  ,  k 

*  *k  "  X 

,  m  *  4) 

MUM 

*  X^e  = 

Y 

<5  =  Mu*  -  til 

1  1  error 

bound  | | 

0 

0.000  069 

0.008 

3 

0.5 

0.000  114 

0.007 

3 

0.75 

0.000  157 

0.006 

3 

0.875 

0.000  192 

0.005 

4 

TABLE  3  (w  r  1,  m  *=  4, 


<k  as  in  (6,6)) 
k 


c 

6  -  Mu*  -  til 

x  1  1 Z  I  I  =  X^e  = 

;  |  error  bound  1  1 

1 

0.000  077 

0.008  1 

0.5 

0.000  088 

0.007  9 

0.2 

0.000  116 

0.007  3 

0.01 

0.000  256 

0.006  3 
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It  was  not  necessary  to  use  orthogonalized  functions 
in  order  to  avoid  undesired  consequences  of  rounding  errors  in 
solving  the  linear  system.  (For  corresponding  error  estimations 
for  partial  differential  equations,  we  had  to  use  orthogonalized 
functions 

6.3  CONCERNING  THE  ACCURACY  OF  THE  ERROR  ESTIMATION 

For  the  examples  of  Section  6.2,  the  actual  error  |u*  -  <f>| 
is  much  smaller  than  the  error  bound,  and  the  ratio  6/ A  decreases 
with  increasing  m.  This  has  to  be  expected  in  general,  if  one 
calculates  error  bounds  for  different  m  in  the  same  simple  way,  as 
we  did. 

If  one  requires  that  the  defect  is  orthogonal  to  all  polynomials 
up  to  the  degree  m  -  1,  or  to  functions  with  similar  behavior,  the 
defect  assumes  the  value  0  at  least  m  times  in  the  interval  [0,1], 
The  error  u*  —  <p  then  usually  oscillates  similarly.  However,  because 
the  defect  d[0]  *  L[u*  -  0]  contains  derivatives  of  the  error,  this 
defect  will  not  decrease  as  fast  as  the  error,  with  increasing  m.  For 
example,  if  L[u]  =  -u",  and  u*  =  f(m)  sin(m-l)TTx  with  some 

A  fS 

(decreasing)  function  f(m),  then  d[4>m]  =  (m-1)  tt  (u*-4>m) . 

If  one  wants  to  get  an  error  bound  which  is  close  to  the  actual 
error,  this  bound  need  have  a  similar  oscillatory  behavior  as  u*  -  <f>. 

To  calculate  such  a  bound  requires  much  more  work  than  to  get  such  simple 
bounds  as  in  Section  6.2.  This  additional  effort  is  better  used  to 
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obtain  better  approximations. 

For  example,  let  <p  *  be  the  approximations  of  Case  1  for  different 

m  *  3, 4, 5, 6,  and  denote  by  A  the  corresponding  value  A  in  Table  1. 

m 

3 

Then,  for  4>  *  $  ,  the  inequality  (6.3)  is  satisfied  for  A  *  1  and 

z  -  c  with  c  ■  +  A^z°  and  z°  as  in  (6.5).  This  function  ; 

is  a  bound  for  j  u*  -  <f>  |  which  is  not  much  larger  than  |u*  -  <J>  | 

3  6  0 

itself,  and  the  similar  function  <p  -  <p  +  A^z  is  an  even  better  bound. 

3 

But,  it  is  not  reasonable  to  use  the  estimation  |u*  -  $  |  <_  c, .  It  would 
be  better  to  use  |u*  -  <f>^|  <_  A^z°. 

The  question  of  practical  interest  is  not,  how  close  to  the 
actual  error  the  error  bound  is,  but  how  snail  a  bound  one  can  get  with 
a  certain  amount  of  work. 

If  one  wants  to  get  exact  error  bounds,  one  has  to  take  into  account 
the  rounding  errors.  The  rounding  errors  which  occur  during  the 
calculation  of  the  constants  a ^  might  influence  the  accuracy  of  the 
approximation  <P ,  but  they  do  not  destroy  the  validity  of  the  correspond¬ 
ing  error  estimation.  Only  the  rounding  errors  which  occur  when  <t>  is 
calculated  with  the  computed  a^,  and  those  which  occur  during  Step  E 
have  to  be  considered  in  some  way.  In  particular,  one  has  to  be  careful 
in  calculating  the  defect  d  [ ]  which  usually  is  a  small  number  computed 
as  a  difference  of  larger  numbers.  Eventually,  one  has  to  compute  the 
defect  with  double  precision.  This  was  not  necessary  for  our  examples. 
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6.4  EXAMPLES  OF  THE  FOURTH  ORDER 

We  treated  the  boundary  value  problem 
TV 

u  -  bu"  +  P(1  +  x(l  -  x)]u  *  1000(1  +  x)  (0  x  1) 
u(0)  -  u'(0)  -  u(l)  -  u'(l)  -  0 


for  different  constants  b  and  V.  Approximate  solutions  4*  have 
been  calculated  using 


<t>0(x)  =  0,  4>1  -  x2(l  -  x)2xi  1  (i  -  1,2, .. .  ,m) , 

w(x)  =  1,  -  xi  1  (1  -  1,2, 


The  method  of  error  estimation  in  Section 
(5.22)  is  satisfied  for  c(x)  =  y[l  +  x(l-x)], 


6.2  can  be  applied  if 
i.e.  if 


-IOtt  <  b  < 


jA1(b)  <  u  <  jA*(b) 


(6.6) 


If  b  and  y  are  small  enough,  one  can  use 

2  2 
Z  ■  X  (1  -  x) 


(6.7) 


because  then  L[z](x)  >  0.  This  was  done  in 
CASE  1:  b  -  0;  y  -  1. 

The  Table  4  contains  some  of  the  corresponding  results: 


$(4)*  ||d[$]||,  and  A||z||  with  A  *  | |d[$]/L[z] | | . 


(6.8) 
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The  values  $(4)  are  listed  with  all  decimals  that  were  calculated 
by  the  computer.  Therefore,  the  last  decimals  might  be  influenced 
by  rounding  errors. 

For 

CASE  2:  b  *  6,  u  -  100, 

different  functions  z  have  been  used,  namely  z  in  (6.7),  and  the 
solutions  z  =  z^  of  the  problem 


zIV  -  bz"  +  k*pz  *  1,  z(0)  ■  z'(0)  ■  z(l)  *  z'(l)  *  0 


for 


k  -  1,  1.25,  2,  4,  10.  (6.9) 

Except  the  cases  m  =  4,  k  =  8;  m  =  4,  k  =  10;  and  m  =  7,  k  *  10, 
the  maximum  of  d[4>]/L[z]  occurred  at  one  of  the  boundary  points, 
so  that  then 


| error  bound | |  =  A | | z | 


IdUll  1 
L[z](0) 


v|  |d[<|>]  |  | , 


(6.10) 


where  the  numbers  v  are  listed  for  different  cases  in  Table  5A. 


The  Table  5B  contains  the  quantities  (6.8)  for  Case  2  with 


For 


CASE  3: 


b  =  20;  \i  =  200,  400,  800,  1200, 
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we  have  checked  how  the  approximation  and  estimation  procedure  behaves 
with  increasing  c(x).  For  b  *  20,  the  second  relation  in  (6.6) 
becomes 

595  *  yA1(b)  <  v  <_  jA*(b)  «  1229 

Among  the  functions  z^  with  k  as  in  (6.9),  z^  gives  the  smallest 
bounds  for  most  of  the  numbers  m  which  we  used.  Table  6  contains 
some  of  the  corresponding  results. 

Finally,  in 

CASE  4:  b  -  20;  y  «  -200,  -400,  -500, 

we  checked  how  the  error  bounds  behave  with  decreasing  c(x).  Some 
results  are  given  in  Table  7.  For  simplicity,  we  used  z  *  z^ 
although  other  functions  z^  might  give  better  bounds. 
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TABLE  4  (b  -  0,  u  «  1) 


m 

_ Hh) _ 

IldUlll 

2 

3.899 

673 

1 

2.595 

3 

3.896 

778 

0 

1.258 

4 

3.896 

778 

4 

1.232 

5 

3.896 

572 

6 

0.117 

TABLE  5A  (b  »  6,  y  «  100) 


z 

v  in 

(6.10) 

z  as  in  (6.7) 

209 

/  1 

925 

l  1.25 

855 

z.  with  k  =  (  2 

k  \ 

673 

I  4 

325 

\  10 

809 

TABLE  5B  (b  =  6,  p  =  100) 


m 

_ ftps) _ 

IldUlll 

3 

2.797 

274 

9 

104.98 

4 

2.797 

275 

1 

51.27 

5 

2.791 

040 

6 

4.67 

6 

2.791 

041 

5 

2.47 

7 

2.790 

992 

4 

1.32 

X  I  z  |  I  *  || error  bound  | 


TABLE  6  (b  -  20) 


A | | z  |  |  •  ||  error  bound 
for  z  ■  z. 


1200 


1.551 

691 

4 

241.09 

0.13 

1.551 

691 

4 

57.51 

0.032 

1.556 

082 

7 

0.904 

0.000 

67 

1.556 

083 

9 

0.602 

0.000 

66 

1.102 

873 

6 

214.90 

0.083 

1.102 

873 

6 

i 

22.61 

0.071 

1.102 

733 

9 

20.33 

0.064 

1.102 

733 

9 

12.74 

0.049 

1.102 

88  9 

5 

1.40 

0.000 

44 

0.851 

844 

09 

262.68 

0.076 

0.851 

844 

12 

63.02 

0.014 

0.850 

614 

11 

39.51 

0.008 

6 

0.850 

614 

12 

25.55 

0.007 

1 

0.850 

904 

94 

2.54 

0.000 

55 

55 
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